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INTRODUCTION 


I . 

The  interior  ballistic  models  to  date  have  captured  the  overall 
performance  of  the  regenerative  liquid  propellant  gun.  Accordingly,  recent 
interest  has  shifted  toward  a  more  complete  treatment  of  the  regenerative 
liquid  propellant  gun  interior  ballistic  process.  In  an  attempt  to  describe 
one  of  the  details  of  the  process,  we  address  the  distributions  of  velocity 
and  pressure  in  the  liquid  reservoir  of  the  regenerative  liquid  propellant 
gun.  The  present  work  is  a  continuation  of  the  treatment  of  the  liquid 
injection  process  presented  by  the  authors1  *  in  which  a  simplified  treatment 
of  pressure  distribution  was  considered  as  part  of  the  injection  model.  The 
distributions  are  derived  from  a  modified  Lagrange  distribution  with  area 
change  to  account  for  the  shape  of  the  regenerative  piston  and  the  center 
bolt.  The  results  are  then  compared  to  a  one-dimensional  simulation. 

II.  CONTROL  VOLUME 

The  regenerative  liquid  propellant  gun  process,  illustrated  in  Figure  1, 
is  initiated  by  firing  a  primer  which  pressurizes  the  initial  combustion 

chamber.  The  chamber  pressure  acting  or.  the  piston  forces  it  to  the  rear, 
compressing  the  liquid  propellant  in  the  reservoir.  After  an  initial 

transient  period,  the  pressure  m  the  liquid  reservoir  will  be  greater  than 
the  combustion  chamber  pressure  due  to  the  differential  area  across  the 

infection  piston,  resulting  in  the  injection  of  liquid  propellant  into  the 
combustion  chamber. 

In  this  paper  we  are  concerned  only  with  the  liquid  propellant 

reservoir.  The  combustion  chamber  pressure  is  taken  from  an  experimental 
firing  of  a  30-mm  Concept  VI  fixture  (shown  in  Figure  1),  and  input  as  a 

boundary  condition.  The  control  volume  is  shown  in  Figure  2.  The  contours  of 

the  piston  and  the  reservoir  are  approximated  by  straight  line  segments  as 
indicated.  The  center  bolt  and  the  transducer  block  are  fixed  in  the 

reference  frame  of  the  chamber.  The  origin  of  the  coordinate  system,  fixed  in 
tie  chamber  frame  of  reference,  is  at  the  rear  (left  hand)  end  of  the 

rtservcir,  and  x  is  the  coordinate  along  the  bolt  as  shown  in  Figure  2.  The 
pi  ton  moves  rearward  with  a  velocity  up(t),  and  the  points  s^t),  s2(t),  and 
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III.  BACKGROUND 


A  simplified  treatment  of  the  control  volume  detailed  above  which  couples 
the  piston  motion  and  the  liquid  velocity  was  described  in  an  earlier 
paper.  The  goal  was  to  describe  the  injection  process  by  coupling  the 

motion  of  the  regenerative  piston  and  the  liquid  injection  thereby  capturing 
the  time-dependent  nature  of  the  discharge  coefficient. 

A  simulation  of  a  firing  of  a  Concept  VI  30-mm  test  fixture  at  BRL  had 
indicated  that  the  model  exhibited  a  mean  value  of  the  discharge  coefficient 
which  agreed  quite  well  with  that  of  the  experimental  data  in  the  steady  state 
operation,  but  there  was  lack  of  agreement  in  the  early  values  and  in  the  slow 
rise  to  steady  state.  Since  the  model  had  not  addressed  the  Belleville 
springs,  which  are  utilizec  in  the  test  fixture  to  permit  the  piston  to  clear 
the  seal  on  the  nose  of  the  center  bolt,  the  springs  were  included  in  the 
model  with  results  shown  in  figures  3,  4  and  5.  The  comparison  between 

experimental  and  predicted  piston  motion,  with  the  inclusion  of  the  Belleville 

springs,  is  shown  m  figure  3.  The  actual  piston  travels  about  0.55  cm, 

hesitates  as  the  rear  transducer  block  comes  to  a  stop  against  the  Belleville 
springs,  and  then  accelerates  and  smoothly  completes  its  stroke.  The 
simulation  shows  the  same  qualitative  motion  as  the  actual  piston.  However, 
the  simulated  piston  position  passes  through  the  transient  stage  more  quickly, 
with  the  hesitation  occurring  approximately  1  ms  before  that  in  the 
experiment.  The  reason  for  the  discrepancy  is  not  clear.  In  figure  4  a 

similar  result  is  reflected  m  the  liquid  pressure.  Qualitatively  the 

simulation  reflects  the  early  oscillations  in  experimental  liquid  pressure, 
and,  after  an  initial  transient  period,  remains  in  close  agreement.  However, 
the  oscillations  occur  too  early  in  the  simulation.  The  comparison  between 
the  discharge  coefficient  derived  from  experimental  data  and  that  predicted  by 
the  model  is  shown  in  figure  5.  Although  the  mean  values  at  steady  state  are 
in  reasonable  agreement,  the  model  does  not  capture  the  slow  rise  to  steady 
state . 
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Figure  5.  Comparison  cf  Discharge  Coefficient  vs  Time  from  Simulation 
(Petted  Line!  anc  Experimentally  Derived  (Solid  Line) 


In  order  tc  exanne  the  discrepancies  between  the  experiment  and  the 
mode*,  several  simplifications  ir.  the  model  were  examined  and  removed.  The 
theoretical  model  ear.  be  broken  into  two  ma;or  components.  The  first  involves 
the  derivation  cf  equations  describing  the  distributions  of  velocity  and 
pressure  m  the  liquid  reservoir  from  a  continuum  analysis.  The  second 
involves  the  coupling  of  these  equations  into  the  momentum,  equation  for  the 
piston.  In  this  paper  the  first  stage  of  the  derivation  is  discussed,  namely 
the  distributions  of  velocity  and  pressure  in  the  liquid  reservoir.  In  the 
earlier  model  the  contours  of  the  outer  piston  and  inner  bolt  were  not 
included  in  the  representation  of  the  pressure  distribution  in  the  liquid 
reservoir;  they  are  now  explicitly  considered.  It  is  not  possible  to  evaluate 
the  impact  of  these  additional  considerations  on  the  comparison  with 
experimental  data  until  the  second  stage  of  the  derivation  is  completed. 
Therefore,  to  assess  the  validity  of  this  stage  of  the  modeling,  a  comparison 
is  made  to  a  one-dimensional  simulation  of  the  liquid  reservoir. 
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IV.  EQUATIONS  OF  MOTION  FOR  THE  FLUID  AND  VELOCITY  DISTRIBUTION 


The  equations  of  motion  for  the  fluid,  continuity  and  momentum  equations, 
are  written  to  include  area  change  as  the  piston  moves  rearward,  we  note  that 
the  area  through  which  the  fluid  flows  is  a  function  of  both  time  and  position 
since  the  contoured  piston  moves  rearward  over  a  contoured  bolt.  The 
equations  of  motion  are  then 

3/  3t  (  pA)  +  3/  3x  (  pvA)  =0  (i) 

3/  3t  (  ova)  +  3/3x  (  pv2A)  =  -  A  3P/3x  .  (2) 

We  obtain  the  equation  for  the  velocity  distribution  from  the  continuity 
equation  (1).  From  (1),  using  the  Lagrange  approximation,  i.e.  assuming  that 
the  density  is  spatially  uniform  at  each  timestep,  3p/ 3x  =  C,  we  have 

A  (1/p)  3p/3t  =  -  3A/3t  -  3/3x  (vA)  .  (3) 

The  goal  is  to  eliminate  the  partial  derivatives  with  respect  to  t,  time,  and 
recast  the  equation  in  3/ 3x  .  The  resulting  equation  is  then  integrated  with 
respect  to  x  to  obtain  an  expression  for  velocity  as  a  function  of  time  and 
position.  Consider 


(4) 


where  mL  is  the  mass  of  liquid  in  the  reservoir,  and  VR  is  the  volume  of  the 
reservoi r . 

Now, 


mL  =  -  CA3V3  -  %A3 


and 


(5) 
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(6) 
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V 


such  that, 


1  if 

p  3t 


[V3A3 


u  A  1 
P  R 


(7) 


where  v3  is  the  velocity  of  the  exiting  liquid,  A3  is  the  exit  vent  area,  AR 
is  the  projected  area  of  the  piston  on  the  liquid  reservoir  side,  and  Up  is 
the  velocity  of  the  piston.  The  exit  vent  area,  A3,  is  the  annular  area 
between  the  piston  at  S3  and  the  bolt.  A3  is  time  dependent,  but  it  is  not  a 
function  of  position  x.  The  time  rate  of  change  of  mass  includes  both  the 
mass  flux  through  the  orifice  with  respect  to  the  bolt  and  the  additional  mass 
flux  intc  the  chamber  due  to  the  piston  motion.  It  is  notea  that  we  have 
chosen  a  laboratory  reference  frame  stationary  on  the  bolt  with  0.0  at  the 
back  wall  and  the  positive  x-axis  forward.  Since  the  natural  motion  of  the 
liquid  is  to  the  right  through  the  ver.t,  and  the  natural  motion  of  the  piston 
is  tc  tne  left  toward  the  back  wall,  v3  is  positive  to  the  right  while  up  is 
positive  tc  the  left. 


It  car.  be  shown  that  the  tire  rate  of  change  of  area  is 

'  2 

A!  x ,  t )  =  u^_  3  3x  (  T?  (  x ,  t ) ) 

where  R(x,t)  is  the  radius  cf  the  inner  contour  of  the  piston  where  we  have 
explicit}-  noted  that  area  is  dependent  on  position  and  time. 


Equation  (3)  then  becomes 


-r-  <va)  = 

ox 


3[  1TR  (  x ,  t)  ] 
Up  3x 


v  A  -  u  A 

L.U....E.5  )  «Ilt, 

VR 


(8) 


Integrating  with  respect  to  x  and  noting  that  v(0,t)  =  0  we  obtain 


v(  X,  t)  =  u 


2  2 

t;f  -  R  t  x  ,  t ) 
A?x, t) 


•  V3/_;W  Y(x.t) 

v”  A?x7tT 


(9) 
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where  V(x,t)  is  the  volume  of  the  reservoir  from  the  back  wall  to  the  position 
x.  Equation  9  defines  the  velocity  distribution  in  the  liquid  reservoir. 

It  is  useful  to  rewrite  equation  (9)  in  terms  of  the  liquid  areas.  Note  that 

it[R2q  -  R2  ( x ,  t)  )  *  *[(R20  *  “  (R2(x,t)  -  rfa2(x))]  *  \  ~  A(x,t) 

where  we  have  used  the  fact  that  the  radius  of  the  bolt  varies  by  only  a  small 

amount  over  its  length,  e.g.  the  maximum  radius  is  1.79  cm  and  the  mimimum 
2  2 

1.65  cm.  Thus  r  -  r  (x)  is  a  very  small  number.  Using  this  approximation, 
0  b 

equation  (9)  becomes 


u  [A,.  -  A(  x ,  t )  ]  ( v  A  -  u  A  ) 

p  W  3_ 3  __p  R_  Vtx^t^ 

"Z  ~  ATxTtT  +  V~  A(x7t” 


V.  PRESSURE  DISTRIBUTION 


Using  equation  (1)  in  (2)  and  integrating  with  respect  to  x  from  O  to  x 


we  obtain 


P ( x , t )  =  P  (t)  -  1 / 2  c  v  ( x , t )  -  p  ,  *  v(x’,t)  dx'  (11) 

O 

in  which  p^(t)  is  the  breech  pressure  at  a  given  time.  Integrating  the  time 
derivative  of  (10) 
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The  relationship  of  breech  pressure  to  P(x,t)  is  given  by 


P(x,t)  =  PQ(t)  -  1/2pv  ( x ,  t) 


Combining  (10),  (12),  and  ( 1 9 ) 


p  /  v ( x ' ,t)  dx’  . 
0 
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11 


(21  ) 


pz.,u  s3 

-  pz  u  x  +  -• - -  f  xA(x,t)dx  . 

1  p  V 

R  0 

The  integrals  involve  A(x,t),  the  area  at  x  at  a  given  time,  ard  V(x,t), 

the  volume  of  the  reservoir  from  the  back  wall  to  x  at  a  given  time.  The 

integration  is  performed  numerically. 

In  order  to  numerically  integrate,  we  first  develop  expressions  f  ir  the  area, 
A( x , t ) ,  and  volume,  V(x,t). 

Approximating  the  contour  on  the  inner  surface  of  the  piston  by  the 
dashed  straight  line  as  indicated  in  Figure  1,  we  can  express  the  radius  of 
the  piston  as  measured  from  the  center  line  as 

V  *1 

R(x,t)  =  {  R1  +  ------  (x  -  s^)  [  1  -  H(s1  -  x)]} 

H(S2  -  x)  +  R  l  1  -  K(s2  -  x)]  H(S;J  -  >i  (22) 

where  x  indicates  the  position  or.  the  bolt,  R1  indicate?  the  radius  of  the 

piston  at  s1  and  R2  is  the  radius  of  the  piston  at  s2.  Kx)  is  the  Heaviside 
function  defined  as  H(x)=0  for  x<0  and  H(x)=i  for  x>  =  0  The  radius  of  the 
inner  surface  of  the  piston  at  a  position  x  is  time  deper  dent  since  the  piston 
is  moving. 

Similarly,  the  radius  of  the  bolt  at  any  position  x  on  the  bolt  can  be 
expressed 


r  (x) 
D 


(X  -  x  )  [  1 


H(x 


x) ]  }  H(x2  -  x)  + 


r^  (  1  -  H(x2  -  x ) ]  H(x3  -  x)  . 
The  cross-sectional  area  of  the  liquid  is  then  giver,  by 


(23) 
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(24) 


2  2 

A( x, t)  =  *R  (x,t)  -  irr,  (x) 

D 

and  the  volume  of  liquid,  V(x,t),  is 

V(x,t)  *  /*  A( x ' , t)  dx’  .  (25) 

0 

Returning  to  (13),  although  it  is  possible  to  obtain  an  analytic  solution 
for  some  of  the  integrals,  the  representation  is  dependent  upon  the  location 
in  the  reservoir  and  the  particular  geometry  of  the  problem.  Thus,  to 
maintain  the  generality  of  the  solution,  it  was  decided  to  obtain  the  values 
of  the  integrals  numerically.  This  procedure  yields  satisfactory  results,  and 
it  was  felt  that  no  further  refinement  was  necessary.  The  pressure 
distribution  is  then  given  by  (21)  with  the  velocity  v(x,t)  given  by  (10). 

VI.  RESULTS 

The  velocity  and  pressure  distributions  predicted  by  the  code  are 
reasonable  fror  a  flui.d  dynamics  viewpoint,  namely,  a  constriction  causes  the 
velocity  to  rise  and  the  pressure  to  fall.  In  particular,  the  calculated 
distributions  reflect  the  geometry  of  the  30-mm  test  fixture  at  the  Ballistic 
Research  Laboratory. 

The  predicted  velocity  and  pressure  distributions  can  be  classified  into 
one  of  two  cases.  For  the  majority  of  the  injection  cycle,  the  piston  head  is 
over  the  flat  portion  cf  the  bolt  as  shown  in  Figure  6  with  points  A,  B  and  C 
denoting  fixed  locations  on  the  piston  as  shown.  The  slant  section  of  the 
piston  is  between  A  and  B  with  the  flat  vent  from  B  to  C.  A  typical  pressure 
distribution  computed  by  the  model  is  shown  in  Figure  7  at  3.0  ms  where  the 
piston  head  is  positioned  over  the  flat  section  of  the  bolt.  The  pressure 
distribution  follows  the  shape  of  the  piston  as  expected;  relatively  constant 
along  the  flat  portion  of  the  piston,  decreasing  over  the  converging  slant 
section  of  the  piston  from  A  to  B,  and  leveling  out  again  through  the  vent 
from  B  to  C.  The  slight  increase  in  pressure  of  approximately  1  MPa  may 
indicate  difficulty  with  the  numerical  integration  or  differentiation  at  some 
timesteps.  From  the  Bernoulli  equation,  we  would  expect  the  velocity 
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distribution  to  generally  follow  the  pressure  distribution  inversely,  i.e.  as 
the  pressure  decreases,  the  velocity  increases.  Figure  8  shows  the  velocity 
distribution  (expressed  in  scientific  notation)  at  3.0  ms  with  the  expected 
result;  the  velocity  gradient  is  shallow  from  the  back  wall  of  the  reservoir 
to  the  beginning  of  the  slanc  section  at  A,  increases  from  A  to  B  on  the 
converging  slant  section  of  the  piston,  and  maintains  a  relatively  constant 
velocity  through  the  vent.  Although  the  velocity  appears  to  be  zero  over  the 
flat  section  of  the  piston,  a  more  detailed  look  at  the  velocity  distribution 
shows  a  linear  increase  from  zero  at  the  back  wall  to  approximately  60  cm/s  at 
position  A,  the  beginning  of  the  slant  section  of  the  piston. 


Figure  7.  predicted  Pressure  Distribution  in  the  Liquid  Reservoir  at  3.0  ms 
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A  second  case  of  piston  position  relative  to  the  bolt  is  shown  in  Figure 
9.  During  the  startup  cycle  the  piston  is  positioned  such  that  the  vent  from 
B  to  C  passes  over  the  slant  section  of  the  bolt  from  D  to  E.  Thus,  exiting  I 

liquid  encounters  a  converging  area  from  A  to  B,  a  flat  section  of  the  vent, 
and  an  additional  converging  area  when  C  is  between  D  and  E.  It  is  expected 
that  the  predicted  pressure  and  velocity  distributions  will  reflect  the 
additional  area  change.  The  pressure  distribution  at  1.3  ms  is  shown  in  | 

Figure  10.  Again  the  pressure  is  stable  over  the  flat  section  of  the  piston, 
drops  between  points  A  and  B,  and  flattens  through  the  vent.  However,  as  the 
liquid  encounters  the  convergent  area  on  the  bolt  from  D  to  E,  the  pressure 
distribution  shows  an  additional  drop.  The  velocity  distribution  in  the  | 

reservoir  at  1.3  ms,  Figure  11,  again  shows  the  expected  relationship  between 
the  velocity  and  pressure  distributions. 

The  magnitude  of  the  velocity  and  pressure  gradients  at  1.3  ms  is  lower  I 

than  at  3.0  ms,  reflecting  the  lower  piston  velocity  near  the  beginning  of  the 
injection  cycle. 


I 


li 


Figure  9.  Relative  Position  of  the  Piston  and  the  Bolt  for  the  Beginning 

Portion  of  the  Injection  Cycle 
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(MPa> 


Figure  i C .  Predicted  Pressure  Distribution  in  the  Liquid  Reservoir  at  1.3  ms 


VII.  COMPARISON  TO  ONE-DIMENSIONAL  MODEL 


A  comparison  of  the  predicted  velocity  and  pressure  distributions  to  a 
one-dimensional  model  developed  by  Coffee^  gives  a  quantitative  assessment  of 
the  model.  The  one-dimensional  representation  differs  from  the  model 
described  in  this  paper  in  that  the  back  wall  is  moved  to  simulate  piston 
motion  as  taken  from  an  experimental  firing  of  the  30-mm  test  fixture  at  the 
BRL.  Since  the  velocity  of  the  back  wall  is  non-zero  in  the  one-dimensional 
model,  the  piston  velocity  was  subtracted  from  the  1-D  velocity  distribution 
for  purposes  of  direct  comparison.  The  piston  positions  were  matched  with 
respect  to  the  fixed  coordinate  system  on  the  bolt  at  a  time  of  3.4  ms  in  both 
models . 

The  comparison  of  the  velocity  distributions  is  shown  in  Figure  12. 
Quantitatively  the  two  models  show  good  agreement  on  the  flat  section  of  the 
piston  and  the  vent  with  a  difference  in  exit  velocities  of  about  1%.  The 
major  discrepancy  between  the  models  occurs  on  the  slant  section  of  the 
piston.  This  is  felt  to  reflect  the  approximation  of  the  partial  differential 
equations  with  ordinary  differential  equations. 

The  comparison  cf  pressure  distributions  between  the  model  and  one- 
dimensional  simulation  is  shown  in  Figure  13.  The  gradients  show  good 
qualitative  and  quantitative  behavior  with  a  difference  of  .2%  at  the  back 
wall. 


Future  work  will  focus  on 
complete  injector  model. 


integration  of 


these  equations  into  the 
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Velocity  (cm/a) 


1.00  3.00  5.  00  7.00  B.  00 


Poa  1 1 1  on  (cm) 


9 


VIII.  CONCLUSIONS 


A  simplification  of  the  continuity  and  momentum  equations  yields 
equations  to  describe  the  velocity  gradient  and  pressure  gradient  in  the 
liquid  reservoir  of  a  regenerative  liquid  propellant  gun.  The  resulting 
equations  are  applicable  to  other  time-dependent  geometries  provided  the 
assumption  that  3p/3x  =  0  is  acceptable.  When  applied  to  the  30-mm  test 
fixture  at  the  BRL,  the  velocity  and  pressure  distributions  are  reasonable 
from  a  fluid  dynamics  perspective  and  compare  favorably  with  results  from  a 
one-dimensional  code.  Future  work  will  concentrate  on  fully  exploiting  the 
pressure  distribution  equation  in  modeling  the  relationship  between  piston 
movement  and  liquid  flow  from  the  reservoir. 
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LIST  OF  SYMBOLS 


*  A(x, t) 


rb(xl 
R(  x ,  t ) 


2 

Cross-sectional  area  of  the  vent,  cm 

2 

Cross-sectional  area  of  the  liquid,  cm 

2 

Cross-sectional  are  of  the  piston  on  reservoir  side,  cm 

2 

Cross-sectional  area  of  the  piston  at  x,  cm 

2 

Cross-sectional  area  of  the  piston  at  the  back  wall,  cm 

Mass  flux  of  liquid  into  the  combustion  chamber,  gm/s 


Combustion  chamber  pressure,  MPa 

Space  mean  pressure,  MPa 

Radius  of  the  bolt  at  tne  back  wall,  cm 


Radius  of  the  bolt  at  x,  cm 


Radius  of  the  pistor  at  x,  cm 

Velocity  of  the  piston,  cm/s 

2 

Acceleration  of  the  piston,  cm/s 
Velocity  of  the  _iquid  at  S3,  cm/s 

2 

Acceleration  of  the  liquid  at  S3,  cm/s 
Volume  of  the  liquid  reservoir,  cm^ 

Density  of  tie  liquid  at  a  given  time,  gtr./cm 
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Concept  6,  30-mm  fixture- -Rd  8  after  springs 


COMBUSTION  CHAMBER  AREA 
PISTON  AREA- -C  CH  SIDE 
PISTON  AREA- -RES  SIDE 
LENGTH  L  PRIME 
LENGTH  OF  VENT 
PISTON  OFFSET 
PISTON  MASS 
VOLUME  LIQUID 
VENT  OPTION 

STRAIGHT  LENGTH  OF  PIST 
MAX  PISTON  TRAVEL 
DENSITY  LIQUID 
K1 
K2 


44.84700 

34.32600 

23.27800 

1.43200 

1.04000 

0.54400 

2109.20000 

172.62764 

2 

5.94680 

7.37880 

1.43700 

5350.00000 

9.11000 


FRICTION  LOSS  LIQ  OPTION-  0 

FRICTION  LOSS  PIS  OPTION-  0 

TIME-C  CH  PRES  DATA  FILE:  ptoffb4.dat 
GEOMETRY  DATA  FILE:  r8geo55.dat 

GRAPH  DATA  FILE:  30mmc6.gra 

INITIAL  PR  IN  RESERVOIR  -  29.00000 

INITIAL  VEL  IN  VENT  -  0.00000 

INITIAL  PISTON  VELOCITY  -  358.00000 

INITIAL  PISTON  POSITION  -  0.00000 

INTEGRATOR- -TINC  -  0.00010 

INTEGRATOR- -EPS  -  0.00001 

INTEGRATOR- -METH  -  1 

INTEGRATOR- -MITER  -  0 

INTEGRATOR- -KWRITE  -  0 


DIFFERENTIAL  EQUATION  SET:  1 


RAD  PIST3  -  1.83000  RAD  PIST2  -  1.83000  RAD  PIST1  -  3.28000 

RAD  BOLT 3  -  1.79770  RAD  B0LT2  -  1.79770  RAD  BOLT1  -  1.65000 


VOL  FUEL12  -  17.90837 

VOL  FUEL23  -  2.02661 


FLAT  LEN  BOLT  -  0.55880  SLANT  LEN  BOLT  -  0.76000 
I WRITE  -  0 


PR  DISTRIBUTION  OPTION  -  1 

NO.  OF  PTS  FOR  DIST  -  20  IP  -  1 

VEL,  PRES  DIST  DATA  FILE:  distvp.gra 
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